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1 .  INTRODUCTION 


There  are  many  methods  available  to  monitor  atmospheric 
aerosols  near  the  ground  and  in  the  stratosphere  where  vol¬ 
canic  dust  may  remain  for  years.  The  technique  of  interest 
to  this  report  involves  the  measurement  of  neutral  points, 
directions  in  the  solar  vertical  where  skylight  is  unpolar¬ 
ized.  Many  measurements  of  sky  polarization,  especially  the 
position  of  the  neutral  points,  have  been  made  since  1886 
after  the  eruption  of  Krakatoa.  Regarded  as  one  of  the  first 
objective  methods  used  to  study  the  properties  of  atmo¬ 
spheric  aerosols,  the  position  of  the  neutral  points  shows 
seasonal  variations  and  strong  deviations  during  highly  vol¬ 
canic  periods.  Recently,  extreme  results  were  obtained  in 
late  1982  in  Hawaii  following  El  Chichon's  eruption. 

Using  customary  visual  observations ,  Volz  has  monitored 
the  position  of  neutral  points  known  as  the  Arago  and  Babi- 
net  points  from  Bedford  and  Lexington,  Massachusetts .  A 
third  neutral  point  called  the  Brewster  point  is  more  diffi¬ 
cult  to  observe  and  will  not  be  discussed  in  this  report. 
The  data  encompass  an  18  year  period  beginning  in  1968  and, 
therefore,  provide  an  excellent  opportunity  to  study  how 
neutral  points  respond  to  varying  tropospheric  turbidity  and 
layers  of  stratospheric  dust.  Although  not  the  aim  of  this 
report,  these  new  observations  might,  on  the  basis  of  good 
knowledge  of  aerosol  parameters  and  theoretical  tools,  be 


1 


used  to  interpret  the  historical  record  of  neutral  points 
Sekera1,  Jensen2. 

1.1  Organization  of  the  Report 

The  present  effort  is  a  partial  evaluation  of  the  se¬ 
ries  of  observations  made  by  Volz.  Chapter  2  presents  a 
brief  discussion  of  the  nature  of  Arago  and  Babinet  points. 
Chapter  3  presents  the  contents  and  organization  of  the  neu¬ 
tral  point  database — key  material  for  the  present  study. 
Chapter  4  discusses  che  effects  of  tropospheric  turbidity, 
sea  surface  reflections  and  snow  cover  on  the  neutral  point 
data.  Chapter  5  describes  the  method  of  removing  the  effects 
of  tropospheric  turbidity  on  the  position  of  neutral  points. 
These  turbidity  corrections  are  necessary  because  the  move¬ 
ment  of  neutral  points  resulting  from  day-to-day  changes  in 
turbidity  may  overshadow  any  movement  caused  by  strato¬ 
spheric  aerosols.  Until  now,  such  corrections  could  not  be 
made  because  turbidity  was  not  an  accurately  measured  quan¬ 
tity.  In  Chapter  6,  the  neutral  point  data  is  presented  in 
a  manner  which  reveals  their  complex  behavior  in  response  to 
the  injection  of  volcanic  dust  into  the  stratosphere.  It 
specifically  focuses  on  the  eruptions  of  Fuego  in  1974  and 
1978,  and  of  El  Chichon  in  1982.  Chapter  6  also  includes  a 
brief  comparison  of  the  measurements  made  by  Volz  with  those 

1.  Sekera ,  Z.  (1956)  Recent  developments  in  the  study  of 

the  polarization  of  skylight,  in  Advances  in  Geo¬ 
physics  ,  Academic  Press,  Inc.,  New  York,  3:43-104. 

2.  Jensen,  C.  (1957)  Die  Polarisation  des  Himmelslichtes, 

in  Handbuch  der  Geophysik,  Borntraeger,  Berlin, 


8 : 527-620 . 


from  othi*  researchers.  Finally,  Chapter  7  summarizes  our 
results . 

2.  THE  ARAGO  AND  BAB 1 NET  POINTS 

2.1  Explanation  of  the  Neutral  Points 

Arago  and  Babinet  points  result  from  the  interplay  be¬ 
tween  primary  and  higher-order  scattered  sunlight.  Insight 
can  be  gained  into  why  these  neutral  points  are  observed 
through  the  following  argument.  First,  consider  a  single 
scattering  atmosphere  composed  entirely  of  small  spherical 
particles.  Here,  one  would  observe  nearly  unpolanzeci  light 
when  looking  very  close  to  the  sun.  This  is  because  sunlight 
is  unpolarized  and  as  a  result,  the  light  scattered  m  the 
forward  direction  by  spherical  particles  will  be  unpolarized 
also.  Now  consider  an  atmosphere  where  higher-order  scatter¬ 
ing  is  important.  This  time,  an  observer  looking  towards  the 
sun  will  see  light  which  has  been  scattered  at  least  once 
from  angles  other  than  forward  direction.  Because  this  light 
is,  in  general,  partially  polarized,  the  observer  will  mea¬ 
sure  partially  polarized  light  when  looking  near  the  sun. 
Extending  the  argument  a  little  further,  the  addition  of  the 
multiple-scattered  to  primary-scattered  light  (having  a  dif¬ 
ferent  state  of  polarization),  will  lead  to  zero  polariza¬ 
tion  being  observed  at  other  locations  in  the  sky.  Speci¬ 
fically,  the  Babinet  point  is  usually  found  about  20  degrees 
above  the  sun  and  the  Arago  point  (only  observable  when  the 
sun  is  low  in  the  sky),  is  about  20  degrees  above  the  anti- 


solar  point.  Calculations  by  Chandrasekhar  ,  which  incorpo¬ 
rate  higher-order  scattering  into  a  Rayleigh  atmosphere , 
confirm  the  presence  of  Arago  and  Babinet  points,  extended 
calculations  have  been  performed  by  Coulson  et  ul4. 

2.2  Factors  Influencing  Neutral  Points 

I ntu i t 1 vel y ,  one  would  expect  the  Arago  and  Habinet 
points  to  depend  on  the  vertical  distribution  and  the  opti¬ 
cal  depth  of  the  scattering  aerosol.  The  first  factor  will 
be  important  especially  when  a  volcanic  dust  layer  is  [ire- 
sent.  For  a  Rayleigh  atmosphere,  the  effect  of  optical  depth 
on  neutral  points  has  been  investigated  by  Holzworth  arm 
Rao-’  using  the  calculations  by  Coulson  et  a  1  4 .  They  have 
shown  that  both  the  Arago  and  Babinet  points  increase  as  op 
tical  depth  increases.  Consequently,  one  would  expect  the 
position  of  neutral  points  to  depend  on  the  solar  elevation 
angle  because  the  optical  depth  of  the  atmosphere  js  a  t unc¬ 
tion  of  solar  elevation  angle.  In  addition,  assuming  that 
aerosols  have  similar  effects  as  air  molecules,  one  may  ex¬ 
pect  a  seasonal  variation  in  the  position  of  neutral  points. 
This  is  because  a  summer  sky  is  usually  hazier  than  the  win¬ 
ter  sky. 

" 3 .  Chandrasekhar ,  S.  (1950)  Radiative  Transfer,  Oxford 
University  Press. 

4.  Coulson,  K.L.,  Dave  J.V.,  and  Sekera  Z.  (I960)  Tables 

Related  to  Radiation  Emerging  from  a  Planetary 
Atmosphere  with  Rayleigh  Scattering ,  University  of 
California  Press,  Berkley  arid  I.os  Angeles. 

5.  Holzworth,  G.,  and  Rao  C.R.N.  (1965)  Studies  of  sky¬ 

light  polarization,  J .  Opt .  Soc .  Am .  55:403-406. 


Because  nult  j pie  scattering  is  very  much  influenced  by 
reflection  phenomena,  the  presence  of  clouds  and  the  type  of 
ground  suilcice  (snow  or  water,  in  particular)  may  also  af¬ 
fect  t  tie  position  of  neutral  points.  These  effects  cause 
noise  it,  t  no  uata  and  put  restraints  on  interpretation. 

.  i  Ht  asurement  of  Neutral  Points 

;  la-  ut  r  a .  po.it  i.ai.i  were  t  uKcii  at  twc  locations  be- 

uiuii  r  ’  ,n<  >•:  oay  t  <  i  twill  gnr.  and  dawn  vary  greatly 

.  .•  i  ’  r.f  ,  •  a  i  .  ■;  !;.  A.r  h  .  :  ci  Geophysics  Laboratory  (AFGL), 

.  w.e  !  i  t  la-  horizon  were  possible,  was 

•  •  r  r>  •  :  *  w>  f<'!.i.  -  M.f  i  Vd  t  i  ons  during  the  winter, 

v.  i.  i  i  a  .  m  na.ii  t  t.e  Vo  1  z  home.  At  this  loca- 

*.  :  :.i  w>  '  i  ,  I  o .  r.t  s  normally  could  not  be  observed 

:  t  t.e  ;  i  1  f  v .  i  ♦  . .  i.  lt.ql*  v.  a sr.aller  tt.an  7  to  14  degrees. 

.  i  |  -  i  *,):.*  i  .  *  uat  .  ot.i.  ,  suet.  as  t  tie  presence  of  volcanic 

e.eud.s,  a’ii'i'p’  :  wit  ton:-  n,  find  )>etter  locations  for  the 

1  I  r  a  *  .  •  r  ■  . 

"in  :,i  eial  j-,  n.t  s  we  i  «•  found  with  a  visual  pol  ari  scope 
t  t. at  ns*  i  a  tav  ait  [  i  i  si!  as  a  polarizer  and  a  colorless  Po¬ 
laroid  !  i  1'fi  a:,  at.  analyzer.  When  viewed  through  this  in- 

*•  t  rumen t  ,  the  neutral  points  appear  as  a  gap  in  the  polar¬ 
ization  fringes  conimonly  known  as  a  "bridge."  A  pendulum 
witt.  an  elevation  angle  scale  has  been  mounted  on  the  in¬ 
strument  so  that  wiiuii  a  neutral  point  is  in  the  line  of 
Sight,  its  elevation  can  be  read  instantly.  Errors  associ¬ 


ated  with  this  measurement  technique  are  0.5  degrees  for  the 


Arago  point  and  1.0  degrees  for  the  Babinet  point,  but  very 
weak  fringes  sometimes  caused  larger  errors. 


The  primary  data  were  obtained  by  measuring  the  angle 
between  the  horizontal  direction  and  each  neutral  point  as  a 
function  of  time,  that  is,  at  odd  solar  elevation  angles. 
Since  it  is  customary  to  express  the  Arago  and  Babinet 
points  as  the  angle  between  zero  polarization  and  the  anti- 
solar  and  solar  direction  respectively,  the  elevation  angles 
of  the  Arago  and  Babinet  points  were  converted  to  solar  dis¬ 
tances  : 


and 


DA  -  Arago  elevation  angle  +  H 


DB  =  Babinet  elevation  angle  -  H, 
where  H  is  the  unrefracted  solar  elevation  angle  associated 
with  a  particular  observation  of  a  Arago  or  Babinet  point 
(see  Figure  1).  Hereafter,  any  reference  to  the  Arago  and 
Babinet  points  will  be  in  terms  of  DA  and  DB  unless  stated 
otherwise.  DA  and  DB  were  then  plotted  against  the  solar 
elevation  angle  for  angles  between  +20  to  -6  degrees  and  the 
values  at  nine  fixed  solar  elevations  were  determined  by 
linear  interpolation.  For  the  purpose  of  analyzing  the 
data,  it  was  helpful  to  also  calculate  the  difference,  DBA  - 
DB-DA.  Figure  2  provides  an  example  of  these  plots  using 
data  representative  of  nonvolcanic  periods.  The  data  and 
plots  were  then  carefully  checked  for  error.  Since  the  in¬ 
terpolation  was  made  irrespective  of  large  gaps  in  H,  im¬ 
proper  values  were  removed.  Also,  extrapolation  over  not 


H  (DEGREES) 


Figure  2.  Values  of  DA,  DB ,  and 
DBA  at  Fixed  Solar  Elevation  Angles 
for  Selected  Days  Between  January  6, 
1970  and  January  11,  1970.  These 
Data  Are  Representative  of  Nonvolcanic 
Time  Periods.  Note  the  DBA  Scale 
Is  Tv/ice  the  Scale  of  DA  and  DB 


'nw>i»uiTTviii<«u»Wi)itrrWwriiu»ui»iJM»liulMWMwiu(»ji^«jiMWi 


ururw.i 


more  than  2  degrees  in  solar  elevation  was  added  where  war¬ 
ranted  . 

The  new  data  file  for  fixed  H  is  the  subject  of  the 
present  work.  The  choice  of  the  fixed  H  angles  reproduces 
the  course  of  the  original  observations  very  well  except  for 
the  times  when  the  volcanic  influence  caused  DA  and  DB  to 
change  rapidly  as  a  function  of  solar  elevation  angle.  This 
volcanic  effect  on  DA  and  DB  was  observed  only  after  the  ar¬ 
rival  of  the  El  Chichon  dust  cloud  (see  Figure  3).  Near  sun¬ 
rise  and  sunset,  better  agreement  between  the  primary  data 
and  the  interpolated  values  would  have  resulted  if  an  addi¬ 
tional  fixed  solar  angle  of  1.5  degrees  had  been  used.  The 
El  Chichon  effect  also  caused  the  measurements  at  a  given  H 
to  change  considerably  from  set  to  set  until  late  1984. 


3 .  THE  DATABASE 
3.1  Neutral  Points 

The  database  consists  of  approximately  1300  sets  of  ob¬ 
servations  in  all.  Each  set  contains  values  for  DA  and  DB  at 
some  or  all  of  the  following  solar  elevation  angles:  20,  14, 
10,  6,  3,  0,  -2,  -4  and  -6  degrees.  Measurements  were  taken 
both  in  the  morning  (M)  and  in  the  evening  (E),  and  have 
been  noted  accordingly.  The  data  encompass  the  time  period 
extending  from  October  1968  to  June  1986,  but  there  is  a 
large  gap  from  1977  to  1980  when  no  measurements  were  taken. 
Furthermore,  the  observations  are  spurious  in  nature  and 
tend  to  occur  in  groups  of  a  few  clear  days  because  observa- 
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tions  can  only  be  made  in  essentially  cloud-free  conditions. 
Figure  4  shows  the  number  of  morning  and  evening  observa¬ 
tions  as  a  function  of  month  and  year. 

3 . 2  Other  Parameters 

In  addition  to  the  neutral  point  data,  six  parameters 
are  included  which  attempt  to  characterize  the  optical  state 
of  the  atmosphere:  atmospheric  turbidity,  turbidity  type, 
cloud  conditions  and  snow  cover,  solar  aureole  and  the  loga¬ 
rithm  of  the  ratio  of  red  to  green  light  in  the  twilight 
sky.  These  parameters  could  prove  helpful  when  analyzing  the 
neutral  point  data.  Atmospheric  turbidity,  plus  cloud  and 
snow  cover  are  described  below  because  they  are  of  some  in¬ 
terest  to  the  current  work.  The  other  parameters  are  summa¬ 
rized  in  the  Appendix. 

3.2.1  Atmospheric  Turbidity 

In  this  report,  the  term  "atmospheric  turbidity"  will 
refer  to  the  aerosol  optical  density  (B)  at  a  wavelength  of 
500  microns.  The  aerosol  optical  density  is  related  to  the 
aerosol  optical  thickness  by 

B  -  Aerosol  Optical  Thickness/2.30. 

Values  for  B  were  obtained  by  means  of  sun  photometer  mea¬ 
surements  of  the  solar  intensity.  The  measurements  were  gen¬ 
erally  taken  shortly  before  sunset  or  after  sunrise.  In  some 
cases  however,  the  values  had  to  be  estimated  from  noontime 
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Number  of  Observations  Number  of  Observations 


measurements  because  of  cloud  and  horizon  conditions.  For 


values  of  B  less  than  0.1,  the  error  in  the  B  is  0.01  in 
magnitude  and  for  higher  turbidities,  it  is  10%  of  the  value 
of  B . 

3.2.2  Cloud  Condition  and  Snow  Cover 

As  a  rule,  measurements  were  taken  only  when  the  local 
sky  was  nearly  free  of  clouds.  Small  transient  clouds  were 
tolerated  as  long  as  they  were  not  in  immediate  vicinity  of 
the  neutral  points.  The  presence  of  clouds  and  snow  cover 
are  coded  in  the  database  according  to  the  numbering  scheme 
described  in  Table  1 .  A  value  for  the  cloud  condition  was 
chosen  so  that  it  represented  the  average  cloud  cover  ob¬ 
served  during  the  course  of  a  whole  set  of  morning  or 
evening  measurements.  Also,  the  code  numbers  from  50  to  53 
indicate  that  the  snow  cover,  within  10  km  of  the  observa¬ 
tion  site,  was  estimated  to  be  at  least  30%. 

3.3  Organization  of  the  Database 

Before  the  neutral  point  data  could  be  analyzed,  it  was 
necessary  to  update  the  existing  data  with  additional  data 
provided  by  Volz  and  to  then  reorganize  the  database.  The 
first  task  was  done  by  merging  a  data  file  composed  entirely 
of  new  data  with  original  data.  The  new  data  consisted  of 
additional  measurements  of  the  neutral  points  and  the  opti¬ 
cal  parameters  discussed  in  previous  section.  A  unified  data 
format  was  then  developed  because  the  original  data  had  been 
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Table  1.  Cloud  and  Snow  Cover  Codes  Used  in  the  Neutral 
Point  Database 


f — - — -  —  — - - — — 

Code  Number 

1 

* 

Observed  Condition 

00 

No  Clouds 

01 

Clouds  less  than  5  degrees  above  the  horizon 
in  the  direction  opposite  to  the  sun 

02 

Clouds  less  than  10  degrees  above  the  horizon 
in  the  direction  opposite  to  the  sun 

10 

Thin  cirrus  less  than  5  degrees  above  the 
horizon  in  the  direction  of  the  sun 

11 

Thin  cirrus  and  altocumulus 
less  than  10  degrees  above  the 
horizon  in  the  direction  of  the  sun 

12 

More  dense  clouds  less  than  10 
degrees  in  the  direction  of  the  sun 

22 

Clouds  less  than  10  degrees  above  the  i 

horizon  at  all  azimuthal  angles 

33 

Same  as  code  22  plus  scattered  clouds  higher 
up 

34 

Clouds  all  over  sky,  but  less  than  5/8  j 

coverage  1 

40 

1 

Cirrus  and  haze  over  much  of  sky  j 

50 

l 

Snow  cover,  cloudless  skies 

51 

Snow  cover,  plus  the  sky  conditions  in  codes 

11,  12,  and  22 

52 

Snow  cover,  plus  the  sky  conditions  in  code 

33 

53 

Snow  cover,  plus  the  sky  conditions  in  code 

34  j 

*  From  1974  to  1986  no  observations  were  made  during  cloudy 
conditions  described  by  code  numbers  12  to  34,  52,  and  53 


given  in  different  formats.  The  new  format  permits  easier 
access  to  the  data  when  performing  statistical  and  graphing 
tasks.  An  example  of  the  newly  adopted  format  is  given  in 
Table  2.  Each  line  corresponds  to  a  particular  date  when 
neutral  point  data  were  taken.  Missing  neutral  point  data 
are  represented  by  zeros. 

At  this  point,  the  data  were  checked  as  much  as  possi¬ 
ble  for  errors  (i.e.  duplicate  records,  typos,  obviously 
wrong  data,  etc.).  Also,  any  observations  made  at  locations 
other  than  Bedford  or  Lexington,  Massachusetts  were  with¬ 
drawn.  This  is  why  some  line  numbers  do  not  contain  data. 

4.  THE  EFFECTS  OF  VARIOUS  OPTICAL  PARAMETERS  ON  THE  NEUTRAL 
POINT  DATA 

4.1  Atmospheric  Turbidity 

As  mentioned  in  Chapter  2,  the  amount  of  scattering  by 
aerosols  will  affect  the  position  of  the  neutral  points.  To 
determine  the  extent  of  this  effect,  scatter  plots  of  the 
neutral  point  data  as  a  function  of  the  observed  atmospheric 
turbidity  were  constructed.  Separate  plots  were  made  for  DA 
and  DB,  and  for  each  fixed  solar  elevation  angle.  When  mak¬ 
ing  the  scatter  plots,  only  the  data  taken  prior  to  the  ar¬ 
rival  of  the  El  Chichon  dust  cloud  were  used  because  it  was 
clear  that  the  presence  of  the  cloud  altered  the  location  of 
the  neutral  points  much  more  than  changes  in  tropospheric 
turbidity  (and  in  a  way  that  was  different  from  the  turbid¬ 
ity  influence ) . 


Table  2.  Format  of  the  Neutral  Point  Database 
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Figures  5a  and  5b  give  plots  of  DA  and  DB  respectively, 
as  a  function  of  atmospheric  turbidity  for  a  solar  elevation 
angle  of  6  degrees.  For  brevity,  only  the  results  at  this 
solar  angle  are  presented;  similar  trends  in  the  data  were 
observed  at  the  other  solar  elevation  angles(see  Figures  6 
and  7).  Figure  5a  suggests  that  for  turbidity  values  less 
than  0.15,  DA  increases  rapidly  as  turbidity  increases  and 
it  approaches  an  asymptotic  limit  thereafter.  For  turbidi¬ 
ties  greater  than  0.3,  there  is  some  indication  that  DA  de¬ 
creases,  but  there  are  too  few  data  points  to  substantiate 
this  trend.  Figure  5b  indicates  that  DB  also  increases  as  a 
function  of  turbidity,  however,  not  as  much  as  DA. 

4.2  Sea  Surface  Reflection 

For  observations  of  DA  made  at  a  lake  or  sea  shore, 
Fraser^  showed  that  a  decrease  of  up  to  4  degrees  can  be  ex¬ 
plained  if  Fresnel  reflection  from  the  water  surface  is  con¬ 
sidered.  Furthermore,  both  observations  and  theory  agree 
that  such  effects  on  DB  (in  the  morning  observations),  are 
less  likely.  Since  our  data  were  obtained  10  to  20  km  from 
the  Atlantic  Ocean,  one  needed  to  determine  if  the  same  ef¬ 
fects  could  still  be  observed.  Therefore,  morning  and 
evening  observations  of  each  neutral  point  were  compared  be¬ 
cause  the  neutral  point-ocean-sun  geometry  is  (for  most  of 
the  year),  different  at  these  times.  To  avoid  any  confusion 
with  volcanic  effects,  only  the  data  prior  to  1982  was  con- 

b.  Fraser,  R.S.  (1968)  Atmospheric  neutral  points  over 

water,  J.  Opt.  Soc.  Am.  58:1029-1031. 
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sidered.  Furthermore,  those- records  of  data  having  turbidi¬ 
ties  greater  than  U.3  or  cirrus  clouds  present  at  the  time 
of  observation  were  not  included. 

Figures  6a  and  6b  separate  DA  and  DB  according  to 
morning  and  evening  observations  for  a  solar  elevation  angle 
of  14  degrees.  Figure  6a  indicates  that  are  no  systematic 
differences  i  n  DA  lot  the  tv*c  t  i  m«-s  ol  <  i.ser  vut  i  on  .  The  same 
analysis  was  performed  for  the  other  pcs. five  solar  eleva¬ 
tion  angles  and  no  distinct  differences  were  observed.  In 
light  of  this,  it  was  felt  that:  the  Atlantic  was  too  far 
away  to  affect  evening  observations  of  DA  and,  therefore,  no 
corrections  for  it  were  necessary. 

Interestingly,  Figure  6b  suggests  that  cur  morning  ob¬ 
servations  of  DB  may  be  lower  than  those  taken  in  the 
evening.  These  differences  however,  are  too  small  to  be  of 
importance  and  therefore,  morning  and  evening  observations 
of  DB  will  not  be  treated  differently  in  the  remainder  of 
this  report  . 

4 . 3  Snow  Cover 

An  analysis  similar  to  that  in  Section  4.3  was  per¬ 
formed  for  snow-cover  <  d  and  bare-ground  conditions.  This  was 
deemed  necessary  because  the  amount  of  sunlight  reflected  by 
a  snow-covered  surface  is  substantial,  especially  when  com¬ 
pared  with  a  vegetative  one.  Previous  work  hy  Neuberger^ 

7.  Neuberger,  H.  (l')4  1  )  The  i  nt  1  nonce  of  the  snow  cover 
on  the  position  of  Arago's  ncut  n.l  point.  Bull  .  Am . 
Meteorol .  So o.  22:348-  tM. 
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Figure  6.  Morning  and  Evening  Values  of  (a)  DA  and 
(b)  DB  as  a  Function  of  Atmospheric  Turbidity  for 
the  Time  Period  1968  to  1977.  These  Data  Are  for  a 
Solar  Elevation  Angle  of  14  Degrees 
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suggests  that  DA  may  be  lowered  by  about  i  degrees  over 
snow-covered  fields  due  to  the  increased  multiple  scatter¬ 
ing  . 

In  Figures  7a  and  7b,  DA  and  DB  at  a  solar  elevation  of 
20  degrees  are  separated  according  tc  snow  and  non-snow  sur¬ 
face  conditions.  No  systematic  differences  occured  between 
the  two  types  of  surface  conditions.  Similar  conclusions 
have  been  made  by  Chandrasekhar^.  His  calculations  showed 
that  the  position  of  a  neutral  point  is  relatively  insensi¬ 
tive  to  ground  albedo. 

5.  TURBIDITY  CORRECTIONS  APPLIED  TO  THE  NEUTRAL  POINT  DATA 

The  previous  chapter  showed  that  neutral  points  are  not 
influenced  by  the  characteristics  of  the  ground,  only  the 
state  of  the  atmosphere  itself.  That  is,  most  of  the  scatter 
often  observed  in  the  neutral  point  data  can  be  attributed 
to  daily  variations  in  tropospheric  turbidity.  It  seems  rea¬ 
sonable  to  try  to  compensate  or  correct  for  these  day-to-day 
variations  because  they  may  overshadow  any  movement  in  the 
points  caused  by  the  presence  volcanic  aerosols.  This  sec¬ 
tion  describes  the  methods  used  to  compensate  for  tropo¬ 
spheric  effects  during  nonvolcanic  and  volcanic  periods. 
Volz  then  discusses  his  examination  of  the  turbidity-cor¬ 


rected  data. 


Cfcgreer  DA  (in  degrees) 


5.1  Nonvolcanic  Periods 

Using  Figures  5a,  5b  and  similar  plots  for  the  other 
fixed  solar  elevation  angles,  the  first  step  in  the  correc¬ 
tion  scheme  was  to  draw  best-fit  curves  through  the  data, 
and  to  record  the  values  for  DA  and  DB  corresponding  to  a 
turbidity  of  0.05.  Deviations  of  DA  and  DB  from  the  values 
at  a  B  equal  to  0.05  were  then  estimated  for  predetermined 
values  of  B  by  subtracting  the  value  for  DA  and  DB  along  the 
curve  from  the  value  measured  at  0.05.  As  an  example.  Figure 
8  illustrates  how  the  correction,  A  DA,  is  derived  for  an 
atmospheric  turbidity  of  0.2  (see  Tables  3  and  4  for  the 
complete  set  of  corrections  for  DA  and  DB ) .  The  corrections 
were  then  applied  to  the  neutral  point  data  prior  to  the  El 
Chichon  eruption  and  after  May  1983.  This  was  done  by  match¬ 
ing  the  observed  atmospheric  turbidity  for  a  given  set  of 
data  with  the  turbidity  values  in  Tables  3  and  4,  and  by 
applying  the  corresponding  correction  values.  Whenever  an 
observed  turbidity  value  fell  between  two  of  the  turbidities 
in  Tables  3  and  4,  the  correction  was  obtained  by  means  of  a 
linear  interpolation  between  the  corrections  associated  with 
the  two  turbidities. 

With  the  corrections  described  above,  the  neutral  point 
data  are  normalized  in  effect,  to  an  atmosphere  with  0.05 
turbidity.  This  value  of  turbidity  was  chosen  because  the 
majority  of  the  points  are  clustered  around  it.  The  approach 
does  not  take  into  account  other  effects  caused  by  the  other 
parameters  in  the  database  and  therefore,  should  be  treated 


as  a  first  approximation.  Furthermore,  the  scneme  assumes 


that  the  contribution  to  atmospheric  turbidity  from  the 
stratospheric  aerosols  is  negligible,  or  at  least  constant 
with  time.  This  is  probably  true  for  nonvolcanic  periods  but 
not  for  highly  volcanic  periods  such  as  the  months  following 
the  El  Chichon  eruption. 

5.2  Volcanic  Periods 

To  determine  to  what  extent  volcanic  dust  in  the 
stratosphere  affects  the  values  of  atmospheric  turbidity, 
Volz®  provided  the  monthly  averages  of  atmospheric  turbidity 
for  the  period  prior  to  the  El  Chichon  eruption  and  the 
monthly  averages  beginning  in  August  1982  (see  Table  5).  It 
is  clear  from  these  data  that  the  El  Chichon  dust  cloud 
caused  the  atmospheric  turbidity  to  increase  substantially 

during  the  period  from  October  1982  to  May  1983.  Using  mea- 

q 

surements  of  peak  optical  depths,  McCormick  et  al.  have 
monitored  the  northward  movement  of  the  El  Chichon  cloud. 
Their  results  confirm  the  presence  of  the  El  Chichon  dust 
cloud  over  northern  midlatitudes  during  the  same  time  pe~ 
r  i  od  . 


8.  Volz,  F.E.  (1984)  Volcanic  turbidity,  skylight 

scattering  functions,  and  sky  polarizat  ion  ami 
twilight  m  New  England  during  1^83,  Appl ■  Opt.  23: 
2589-2593. 

9.  McCormick,  M.P.,  Fwissler,  T.J.,  Fuller,  W.,  Hunt  H . , 

and  Osborn,  M.T.  (1984)  Airborne  and  ground-based 
lidar  measurements  of  the  El  Chichon  stratospheric 
aerosol  from  90  degrees  N  t«.  56  degrees  S,  Oeof.  Jnt. 
23:187-221. 
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Table  5.  Monthly  Averages  of  Turbidity  at  Lexington 
Massachusetts  for  the  Years  1  980  to  1983. 


Month 

1980  -  1981 

Year 

1  982 

1983 

January 

0 .04 

0.075 

0.11 

February 

0 .06 

0.08 

0.12 

March 

0 . 065 

0 .08 

0.125 

Apr  i  1 

0.075 

0 .09  5 

0 .20 

May 

0.13 

1 

0.13 

0.17 

June 

0.16 

0.17 

0.12 

July 

1 

j  0.18 

1 

0 . 20 

0.11 

August 

1 

0.14 

0.13 

0.10 

September 

0.13 

0.135 

0.11 

October 

0.08 

0.09 

0.05 

November 

0 . 06 

0.10 

0 . 04 

December 

I  0.05 

0.12 

0.05 

Based  on  the  large  increases  of  B,  it  is  inappropriate 
to  apply  the  correction  described  in  Section  5.1  to  the  data 
taken  between  October  1982  and  May  1983.  In  effect,  one 
would  be  overcompensating  for  the  movement  of  the  neutral 
points  because  the  corrections  in  Section  5.1  are  for  tropo¬ 
spheric  aerosols  only.  Therefore,  the  correction  scheme  was 
modified  before  applying  it  to  data  taken  during  the  El  Chi- 
chon  volcanic  period. 

Except  for  one  additional  step,  the  method  used  to  com¬ 
pensate  for  tropospheric  effects  during  the  El  Chichon  vol¬ 
canic  period  was  the  same  as  that  for  nonvolcanic  periods. 
Specifically,  first  the  tropospheric  turbidity  associated 
with  a  given  set  of  observations  was  isolated  by  subtracting 
a  turbidity  value  representative  of  the  El  Chichon  cloud. 
That  is, 

Btrop  “  B  "  V' 

where  Btrop  is  the  tropospheric  turbidity  and  Bv 1  is  the 
turbidity  of  the  volcanic  dust  cloud.  Using  B.  in  place 
of  B,  the  volcanic  data  was  then  corrected  using  the  method 
described  in  Section  5.1  and  the  set  of  corrections  in  Ta¬ 
bles  3  and  4 . 

Monthly  values  for  Bv  1  were  obtained  by  subtracting  the 
average  values  of  turbidity  for  each  month  of  the  El  Chichon 
volcanic  period  ( JTJ )  from  the  average 


values  for  the  same 


months  during  1980  and  1981  (Bnonv).  That  is, 

Bv  Bv  Bnonv 

for  each  month.  Here  it  should  be  noted  that  the  day-to-day 
variations  in  the  El  Chichon  cloud  are  assumed  to  be  rela¬ 
tively  small.  This  seems  reasonable  because  the  cloud  was 
well  mixed  by  the  time  it  reached  the  latitude  of  the  obser¬ 
vations  and  therefore,  any  changes  in  the  turbidity  of  the 
cloud  should  be  slowly  varying  with  time.  The  monthly  values 
of  B  '  are  given  in  Table  6. 


Table  6.  Monthly  Values  of  the  Turbidity  of  the  El 
Chichon  Dust  Cloud  ( •  )  ,  from  Measurements  at  Bedford 
and  Lexington,  Massachusetts 


1 

Oct . 

1982 

Nov . 

Dec . 

Jan . 

Feb . 

1983 

Mar  . 

Apr . 

May 

0.01 

0.04 

0.07 

0.07 

0.06 

0 . 06 

0.12 

0.04 

5.3  Examination  of  the  Corrected  Data 

With  the  average  turbidity  effects  removed,  plots  in 
the  form  of  Figures  2  and  3  were  prepared  for  all  the  neu¬ 
tral  point  data.  Many  probable  defects  were  noted,  espe¬ 
cially  by  looking  at  the  deviations  from  the  seasonal  trend 
of  DBA.  In  most  cases,  Volz  could  determine  the  source  of 
the  error  and  the  data  were  edited  accordingly.  A  number 
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of  errors  also  arose  from  wrong  turbidity  values  and  from 
incorrect  classification  as  not  being  cirrus-veil  cases. 

The  complete  database  has  not  been  included  as  a  part 
of  this  report  because  it  is  too  extensive.  Copies  of  the 
turbidity-corrected  data  can  be  obtained  from  Volz  at  AFGL . 
The  turbidities  reported  with  the  turbidity-corrected  data 
already  have  the  monthly  values  of  Bv 1  subtracted.  If  a  tur¬ 
bidity  equaled  zero  after  subtracting  By <  f  that  turbidity 
has  been  given  a  value  of  0.02. 

6.  PRESENTATION  OF  THE  DATA  BY  GROUP  AVERAGES 

In  Chapter  4  it  was  shown  that  the  neutral  points  were 
not  affected  by  the  characteristics  of  the  ground,  only  by 
the  state  of  the  atmosphere.  As  a  result.  Chapter  5  was  an 
attempt  to  eliminate  the  movement  in  the  neutral  points 
caused  by  tropospheric  aerosols. 

The  corrected  neutral  point  data  is  now  presented  in  a 
manner  that  reveals  their  behavior  during  various  strato¬ 
spheric  phenomena.  Specifically,  the  stratospheric  aerosol 
number  densities  and  observed  twilight  color  ratios  will  be 
used  as  guides  to  separate  the  data  into  groups  that  coin¬ 
cide  with  the  time  periods  of  known  volcanic  influences 
(minor  and  major)  and  seasonal  events.  In  effect,  the  group¬ 
ing  of  sets  of  observations  allows  the  entire  data  set  to  be 
presented  in  a  compact  form  and  still  resolve  any  movement 
in  the  neutral  points  caused  by  a  change  in  the  strato- 


6 . 1  Grouping 

The  grouping  of  the  data  was  based  primarly  on  strato¬ 
spheric  aerosol  data,  but  measurements  of  the  twilight  color 
ratios  { CR )  were  considered  because  they  contain  details 
that  may  help  explain  the  behavior  of  the  neutral  points. 

Figure  9  is  a  contparision  between  the  CR  measurements  by 

Volz  at  Bedford,  Massachusetts  and  the  stratospheric  aerosol 
concentrations  measured  with  a  balloon-borne  dustsonde  by 
Hoffman  and  Rosen^.  The  stratospheric  aerosol  measurements 
were  taken  every  four  to  five  weeks  near  Laramie,  Wyoming 
(42N)  and  have  fluctuations  of  5  to  20%.  Figure  9  indicates 
that  the  CR  measurements  follow  the  same  trends  as  the 

Hoffman  and  Rosen  data,  but  the  CR  are  always  lower  during 
spring  and  summer  than  in  the  fall  and  winter.  The  lower 
values  of  CR  may  be  due  to  the  higher  turbidity  near  the 
tropopause  during  the  spring  and  summer.  Therefore,  in 
addition  to  the  nonvolcanic,  moderate  and  strong  volcanic- 
periods,  groups  of  spring/summer  and  fall/winter  CR 

conditions  were  distinguished  in  the  graphs  and  in  the 
numerical  evaluation.  Details  on  the  40  groups  selected  are 
given  in  Table  7.  Unf ortunately,  the  detailed  trend  of  CR 
during  the  El  Chichon  period  was  not  available  at  the  time 
Table  7  was  prepared  and  therefore,  the  grouping  of  data  in 
the  El  Chichon  period  was  not  optimal. 

The  average  values  and  standard  deviations  of  DA  and 
DB  and  average  values  of  DBA  were  then  calculated  for  each 

10.  Hoffman,  D.,  and  Rosen,  J.  (1986)  Atmospheric  effects, 

SEAN  Bull.  1 1 : No . 5 ,  20. 


Trends  of  Peak  Color  Ratio 


Time  Period 


Group 

Number 


Record  * 

Numbers  Classification 


26 

12/26/74 

2/26/75 

882-898 

WV 

27 

3/3/75 

9/29/75 

899-918 

wv-s 

28 

10/7/75 

1/6/76 

919-926 

(WV) 

29 

1/10/76 

7/25/76 

927-937 

NV-S 

30 

9/6/76 

3/11/77 

938-959 

NV-F 

THREE  YEAR  GAP  IN  DATA 


5/22/80 


32  10/31/81 


5/28/80 


10/2/82  965-995 


33  10/10/82  11/27/82  996-1019 


34  12/4/02  3/26/03  1020-1056 


35  4/23/83  8/25/83  1057-1099  (V-S) 

36  9/1/83  2/27/84  1100-1140  MV 

17  3/1/84  7/25/84  1141-1177  MV-S 

38  8/8/84  1/29/85  1178-1239  (MV) 

39  2/10/85  11/23/8 5  1240-1295  W V-S 

40  12/2/85  6/17/86  1296-1357  ( WV ) 


Bulk  of  Fuego  Aerosol 

Weak  Spring/Summer  Minimum  of  CR 

High  CR 

Decrease  of  CR 
Higher  CR 


Mt.  St.  Helens  Eruptions  On  and 
After  5/18/80;  No  Enhanced  Twilights 
Observed  in  Bedford,  MA 

After  the  Eruption  of  El  Chichon 
from  4/3/82  to  4/4/82  in  Mexico,  the 
Increase  in  Twilight  Effects  Were 
Probably  Due  to  the  Volcanic  Aerosol 
Arriving  in  the  Lower  Stratosphere 

Still  Increasing  Twilight  Effects; 
Measurable  Turbidity  Increases  and 
Significant  Changes  in  Neutral 
Points,  First  in  DB 

Bulk  of  the  El  Chichon  Aerosol;  Some 
Twilight  Displays  Indicate  That  the 
Dust  Layer  Is  Very  High  in  the 
Stratosphere,  Possibly  in  Relation 
to  the  Recovering  of  DB  from  1  to  3 
Degrees  Solar  Depression  Angle. 

Slowly  Declining  Volcanic  Effects 
and  Lower  CR 

Slowly  Declining  Volcanic  Effects 

Further  Declining  Volcanic  Effects 

Last  Traces  of  the  Scattering  by  the 
FI  Chichon  Aerosol  Disappear  from 
t  he  Day t ime  Sky 

Neutral  Points  Nearly  Normal 

Neutral  Toints  Normal 


*  Classification: 


V  -  Strong  Volcanic  Period  and  N**utral  Points  Strongly  Affected 
MV  -  Moderate  Volcanic  Period  and  Neutral  Points  Barely  Affected 
WV  -  Weak  Volcanic  Period  and  No  Effect  on  Neutral  Points 
NV  -  Nonvolcanic  Period 
-S  -  Spr inq/Summer  Minimum  of  CR 
-F  -  Fall /Winter  Maximum  of  CR 
(  )  -  Weak  Effect  s 


of  the  groups  in  Table  7  (see  Tables  8,  y  anil  lU).  When 
computing  the  averages,  data  having  a  turbidity  greater  t nan 
0.30  or  cirrus  clouds  present  at  the  time  of  observat : <  n 
were  not  included.  Furthermore,  some  groups  diu  not  have 
data  at  20  and  -6  degrees  solar  elevation  angle  and  there¬ 
fore,  values  had  to  be  estimated,  'these  values  are  derated 
in  brackets. 

6.2  3-D  Plots  and  Staggered  Plots 

Figures  10,  11  and  12  are  3-D  plots  ut  DA,  PH  and  DBA 
respectively,  as  a  function  of  solar  eievatnr.  ar.gl<  and 
group  number.  The  vertical  axes  have  beer,  chosen  so  that  DA, 
DB  and  DBA  all  have  the  same  scale.  Note  that  the  plotting 

package  lias  assigned  equally  spaced  intervals  between  the 

nine  fixed  solar  elevation  angles  and,  therefore,  rm  ;afa 
at  the  far  left  has  been  compressed.  Alsi  ,  the  perspect.ve 

of  the  3-D  plots  caus*.  s  the  M  Chicle  r.  pa  ak  t<  r,;de  ji.<  s*  of 

the  later  profiles.  By  r«  versing  t  h<  ..equenc*-  of  < ;  r  <  ups  >2 
to  40,  the  slow  decline  of  IB  is  showr.  ir.  Figur*-  i *. 

Staggered  plots  of  da,  [>b  ar,<:  DBA  nr  •  y  r  v  ♦  r.  .  r.  K  ;  y  ..  r  •  s 
14a,  14b  and  14c  for  the  41;  groups.  m<  u  •  n :  u : . .  r  . ,  i. 

of  DA  and  DB  have  been  plot  t ed  us  v«  ;  ♦  :  <  u  1  :  n :  .  :  r  .  y  . .  r  e  s 

14a,  14b  and  14c.  In  contrast  ♦  o  t  h»  i  - 1  flits,  rm 
staggered  pi  ot  s  have  solar  •leva*..;.  •  * .  . ,  *  ;  * 

proper  spacir.g  between  t  r,.  t  :  i  *  !  .  :<  •  . :  .  ,  r  , ,  •, . 

gles  . 
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Table  8.  Averaged  Values  of  DA  and  Their  Standard 
Deviations  (In  Parenthesis)  for  the  Groups  in  Table  7 
The  Brackets  Indicate  Those  Values  That  Have  Been 
Estimated.  Also,  a  Standard  Deviation  of  0.0  Means 
That  Only  One  Observation  Was  Available 


Solai  Klevation  Angle  (Degrees ) 
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Figure  10.  3-D  Plot  of  DA  for  the  Groups  Given  in 

Table  7.  For  Reference,  the  Approximate  Time  Scale  Is 
Given  Along  the  Group  Number  Axis 


Figure  11.  3-D  Plot  of  DB  for  the  Groups  Given  in 

Table  7 


6.3  Discussion  of  the  Plots 


The  3-D  plots  are  more  instructive  than  the  staggered 
plots  because  they  bring  out  general  trends  m  the  data  very 
well.  The  flat  valley  in  DA  shortly  after  sunset  and  before 
sunrise  (i.e.  H  between  0  and  -4  degrees)  is  a  pronounced 
feature  throughout  all  of  the  data.  Prior  to  the  El  Chichon 
period,  the  surfaces  of  DA  and  DB  have  about  the  same 
smoothness,  despite  the  higher  accuracy  of  the  DA  measure¬ 
ments.  Only  the  values  of  DB  in  the  first  group  (if  the 
omitted  observations  1  to  12  in  Table  7  are  included)  seem 
to  have  a  significant  deviation.  This  is  probably  due  to  the 
remnants  of  earlier  volcanic  eruptions. 

In  the  El  Chichon  period,  DB  rises  about  twice  as  much 
as  DA.  The  high  values  of  DB  extend  not  only  to  a  solar  ele¬ 
vation  angle  of  3  degrees,  but  sometimes  to  0  or  -1  degrees. 
Figure  13  indicates  that  the  values  of  DB  were  near  normal 
by  1986,  just  like  the  stratospheric  aerosol  counts  and  CR 
data  ( Figure  9  )  . 

The  surface  of  DBA  is  generally  a  flat  ridge  that 
reaches  a  maximum  at  a  solar  elevation  angle  of  0  degrees. 
During  the  El  Chichon  period  however,  a  peak  protrudes  over 
the  ridge.  Also  the  surface  is  slightly  rougher  than  those 
of  DA  and  DB .  The  edges  are  especially  rough  because  fewer 
observations  were  taken  at  -6  and  20  degrees.  Incidentally , 
the  expression 


DA  +  DB  -40 


produces,  even  for  the  El  Chi  chon  condition,  a  very  similar 
course  as  DBA. 

One  observation  can  be  made  from  the  staggered  plots 
(see  Figure  14):  the  standard  deviations  of  DA  and  DB  m  the 
El  Chichon  period  are  larger  than  those  in  nonvolcanic  peri¬ 
ods.  This  suggests  that  the  day-to-day  movement,  of  neutral 
points  is  greater  during  highly  volcanic  conditions.  No  ex¬ 
planation  for  this  behavior  is  offered  at  the  current  time. 

In  an  attempt  to  find  a  relationship  between  the  de¬ 
tails  provided  by  the  CR  and  the  behavior  of  the  neutral 
points,  a  number  of  pre-El  Chichon  groups  in  Table  7  have 
been  separated  according  to  high  and  low  conditions  of  CR 
(see  Table  11).  Group  a  is  a  comparison  between  neutral 
point  data  taken  after  the  arrival  of  the  Fuego  dust  clouds 
of  1971  and  1974  (high  CR  conditions),  and  background  values 
(low  CR  conditions).  In  groups  b  through  f,  fall/winter  mea¬ 
surements  (high  CR  conditions),  of  the  neutral  points  are 
compared  with  those  for  the  spring/summer  (low  CR  condi¬ 
tions).  Here  it  should  be  noted  that  group  numbers  32 
through  40  in  Table  7  were  not  included  in  the  analysis  be¬ 
cause  these  data  were  dominated  by  presence  of  the  El  Chi¬ 
chon  dust  cloud. 

For  group  a,  no  significant  differences  were  found  be¬ 
tween  the  values  of  DA,  DB  am!  DBA  m  high.  CR  conditions 
and  those  in  low  CR  conditions.  In  light  of  this,  we  t  e<.  1 
that  the  Fuego  dust  clouds  ot  1971  ann  1974  were  probably 
too  diffuse  to  effect  the  position  of  the  neutral  points. 


Supergroups"  for  High  and  Low  CR  Conditions 


In  contrast,  groups  b  through  f  indicate  th.it 


tii*  r  • 


seasonal  movement  of  DA  and  DB .  In  particular,  groups  t, 
and  d  yield  essentially  the  same  differences  betw.a-r. 
fall/winter  values  of  DA,  DB  and  DBA  and  spr ing/ summei  m» a 
surements  (illustrated  for  DBA  in  Figure  15).  Groups  e  and 
f,  on  the  other  hand,  differ  from  groups  b  through  u  in  vui  - 
ious  aspects.  Unfortunately,  we  have  no  explanation  foi  t  la- 
behavior  of  groups  e  and  f  and  consequently,  feel  it  would 
be  premature  to  draw  conclusions  from  the  results  in  Figure 
15. 


Figure  15.  ADBA  as  a  Function  of  Solar 
Elevation  Angle  for  Supergroups  b  through  f 


6.4  Coapanson  with  other  Measurements 

6.4.1  Data  from  Mauna  Loa ,  Hawaii  during  the  El  Chichon 
Period 

'lh*'  h!  -  :  c  in  :■  Oust  cloun  i  «  ached  an  unusually  great 
altitude  over  '  t.«  tropics  shortly  a  ltd  the  eruption.  Lidar 
rrt  urns  [leaked  at  it  to  27  kr> ,  compared  wit.ii  17  to  22  km  on 
most  days  .u  t  t;e  ne  i  tr.fr  r.  m  r  c.  1  a  t  i  t  ude  s  ,  and  the  optical 
thicknesses  w >  ft  up  ti  twice  as  large  disc.  From  his  photo- 
polarmetric  measurements  during  July  and  August  of  1982  at 
Mauna  lea,  Hawa  .  l  ,  Coulson^  has  derived  neutral  point  data. 
H;s  data  exr.i  t  .f  much  more  {-renounced  features  (see  Figure 
;t>;  when  compared  with  the  values  in  Figure  3.  In  particu¬ 
lar,  PH  rose  tt  60  degrees  at  a  solar  elevation  of  4  de¬ 
grees,  descended  momentarily  to  about  0  degrees  at  the  time 
when  the  solar  rays  passed  through  the  aerosol  layer  tangen¬ 
tially  (i.e.  at  a  refracted  solar  elevation  angle  of  0  de¬ 
grees),  and  ascended  again.  To  honor  the  late  Kinsell  L. 
Coulson,  who  first  measured  those  sunset  features,  we  pro¬ 
pose  to  name  them  Coulson  feature's.  In  short,  Coulson  fea¬ 
tures  are-  characterized  by  a  minimum  value  of  DB  near  sunset 
i  or  sunrise)  arid  an  increase  at  negative  solar  elevation 
ar.g  1  es  . 

In  our  observations,  good  examples  of  the  Coulson  fea¬ 
ture  s  seem  to  have  occurred  or. Jy  on  a  few  evenings  between 
octet  ei  .98*.  and  Janiiary  1983.  Four  of  the  best  cases  are 
shown  in  Figure  17;  their  similarity  to  the  Hawaii  results 

i  1  .  Coulson,  K.I..  •  14  8  3)  Effects  of  the  El  Chichon  volcanic 

cloud  ir,  trie  stratosphere  or:  the  polarization  of 
light  from  tin  sky,  Appl  ,  opt  .  2 2 :  1  036-1030. 


DB  (DEGREES)  DA  (DEGREES) 


COULSON  FEATURES 


H  (DEGREES) 


Figure  17.  Examples  of  the  Coulson 
Features  in  Our  Data  for  Selected 
Days  Between  November  6,  1982  and 
December  4,  1982 


is  evident.  Weaker  Coulson  features  were  also  observed  be¬ 
tween  July  and  September  of  1983,  but  under  less  favorable 
observation  conditions.  1 nteresti ng 1 y ,  our  Coulson  features 
seem  to  coincide  with  the  time  periods  when  an  aerosol  layer 
was  observed  above  the  normal  Kl  Chichon  dust  cloud.  The 
aerosol  layer,  which  may  be  the  format  ion  of  sulfuric  acid 
droplets  from  sulfur  dioxide  gas  of  the  Kl  Chichon  eruption, 
has  been  seen  at  altitudes  of  up  to  40  km  (at  northern  mid¬ 
latitudes)  in  the  Kuropean  Jidar  measurements  of  Adriani  et 
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to  compart'  scries  <>t  ohsei  vat  ions  because  lew  are  averaged 
or  are  on  coni  ui  m:  nq  plots.  Therefore,  it  was  desirable  to 
assemble  the  historic  data  r  nt  c  .1  an  .  }  i  <  d  lot  mat. 

In  the  cour  se  t  l  t  in-  ;  s  s<  nt  si  i  h  ,  plot  s  u.  the  lorm  ol 

Figures  2  and  J  were  prcuucea  ‘  <  r  most  of  the  historical 

data.  Depicted  1 1.  the  1  e  i  t  pa  r  t  <1  I  i  g  u ;  e  !  H  are  tlit'  results 

from  Germany  :  dt  uses !  4  ;  t  <  r  the  teats  that  loll  owed  t  he 
eruptions  ol  M  t  .  Sail  rule  ami  M  t  .  h<  :  ee  in  the  s  [  ring  of 
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crease  i  r,  r  lie  value  of  L/B  after  unset  .  Although  not  shown 
in  this  report  ,  1 -’la .  1  sep  h  iitnie;  . :  I  s<  i  w  i  r  e  not  observed  a  f  - 
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7.  Summary 

The  Arayo  and  Baba  net  neutral  points,  defined  as  the 
distances  DA  and  DB  from  the  anti  solar  and  solar  angle  re¬ 
spectively,  were  obtained  at  Bedford  and  I.exington,  Mas¬ 
sachusetts,  for  most  of  the  years  between  1968  and  1986. 
Neutral  point  data  at  nine  fixed  solar  elevations  formed  the 
basis  for  the  present  effort.  Many  observations  extend  to 
solar  elevations  of  20  degrees  and  -6  degrees.  The  data  have 
been  checked  carefully  for  errors  and  combined  with  addi¬ 
tional  measurements  of  atmospheric  turbidity,  turbidity 
type,  cloud  cover,  solar  aureole  and  twilight  color  ratio. 
The  resulting  database  consists  of  a  total  of  about  1300 
morning  and  evening  observation  sets. 

The  main  purpose  of  the  present  effort  was  to  assess 
the  seemingly  undetectable  effects  of  the  two  moderate  vol¬ 
canic  dust  veils  of  Fuego  in  1971  and  1974,  and  the  large 
effect  after  the  El  Chichon  eruption  of  1982.  Prior  to  this, 
the  effect.  of  tropospheric  turbidity  on  the  neutral  points 
was  determined  using  the  data  proceeding  1980.  The  neutral 
points  were  then  corrected  accordingly,  in  addition,  the  ef¬ 
fects  of  snow  covered  ground  and  sea  surface  reflection  of 
sunlight  on  the  neutral  points  were  found  to  be  negligible. 

The  data  set  was  subdivided  into  40  groups  for  plotting 
purposes.  The  grouping  was  based  primarily  on  stratospheric 
aerosol  data,  but  the  details  of  the  twilight  color  ratio 
measurements  were  also  considered.  A  number  of  observations 
were  made  t  1  uni  the  the  dots: 


1)  the  surface  of  DA  contains  a  pronounced  flat  valley 
near  a  solar  elevation  angle  of  -2  degrees 

2)  for  the  pre-El  Chichon  period,  the  3-D  plot  of  DBA 
has  a  much  rougher  surface  than  that  of  DA  and  DB 

3)  the  El  Chichon  effect  on  the  data  took  at  least  three 
years  to  subside 

4)  the  grouping  of  the  data  according  to  high  and  low 
values  of  CR  seems  to  have  yielded  signatures  for  the 
spring  to  fall  effect 

5)  the  weak  eruptions  of  Fuego  were  not  detectable 

6)  the  standard  deviations  of  the  neutral  point  data  in 
the  El  Chichon  period  were  larger  than  those  in  non- 
volcanic  periods 

Our  data  were  then  compared  with  measurements  made  by 
other  researchers.  In  general,  the  increase  of  DB  and  DA 
during  the  El  Chichon  episode  was  smaller  than  it  was  ini¬ 
tially  in  the  tropics.  During  the  Fall  of  1982  and  early 
Winter  of  1983  however,  the  high  values  sometimes  continued 
at  negative  solar  elevation  angles,  with  a  sharp  minimum  in¬ 
tervening  at  actual  sunset  or  sunrise  (that  is,  at  0  degrees 
refracted  solar  elevation  angle).  These  "Coulson  features," 
which  seem  to  have  been  caused  by  the  especially  great 
height  of  the  El  Chichon  dust  cloud  and  the  reduced  volcanic 
turbidity  at  lower  altitudes,  were  not  observed  after  the 
eruptions  of  1902  and  1911,  probably  because  the  high  layers 
of  volcanic  dust  did  not  exist  . 
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Appendix 

Additional  Parameters  in  the  Database 


A . 1  Turbidity  Type 

The  Turbidity  Type  (TT)  is  a  qualitative  observation 
that  gives  a  sensitive  indication  of  the  size  distribution 
in  the  radius  range  from  about  0.2  to  2  microns.  The  obser¬ 
vations  were  made  at  a  solar  elevation  angle  of  about  25  de¬ 
grees  and  require  the  use  of  brown  sunglasses,  both  for  re¬ 
duction  of  glare  and  for  the  supression  of  molecular  scat¬ 
tering.  Cold  air  inflow  is  often  characterized  by  strong 
forward  scattering  and  by  blue  sky  farther  than  three  de¬ 
grees  from  the  sun  (TT=1.3).  By  contrast,  hazier  skies  may 
exhibit  a  nearly  constant  brightness  (and  color)  up  to  20 
degrees  away  from  the  sun  (TT=2.U).  The  turbidity  type  is 
stored  m  the  database  as  TT  x  10. 

A. 2  Solar  Aureole 

Normally,  the  sky  around  the  sun  is  vet >  hi  ignt  due  ti 
the  forward  scattering  by  dust  particles  and  pi  lien.  Them 
particles  typically  have  diameters  between  about  2  and  in 
microns,  comparable  in  size  to  teg  and  clout;  droplets.  As 
shown  by  observations  in  the  course  of  a  clear  day,  the  dust 
and  pollen  are  lifted  by  winti  and  turbulence  1  t  on,  the  sur¬ 
face  industry  and  road  traffic.  They  are  kej  t  at  lint  by 
convection . 


fil 


Rather  crude  estimates  of  aureole  brightness  (A)  were 
obtained  until  1974  by  looking  at  the*  sky  about  0.5  degrees 
from  the  sun,  from  the  shadow  of  a  distant  roof.  If  with  the 
unaided  eye,  the  aureole  was  not  glaring  or  was  barely 
glaring,  the  designation  was  A- 0 .  if  one  sunglass  was  needed 
to  reduce  glare  sufficiently,  then  A=1 ,  etc.  Sky  conditions 
with  cirrus  veils  can  easily  be  recognized  by  this  method 
because  these  ice  crystals  have  very  strong  forward 
scattering  up  to  0.2  degrees  from  the  sun.  The  brownish 
sunglasses  used  to  make  the  measurements  had  an  optical 
density  between  1.0  and  1.3  for  wavelengths  between  0.45  and 
0.55  microns,  and  0.5  at  wavelengths  of  0.64  and  0.35 
microns . 

A . 3  Color  Ratio 

The  logarithm  of  the  color  lutio  i  l.CR )  is  derived  from 
continuous  recordings  duiii.g  twilight  of  the  sky's  intensity 
at  wavelengths  cf  0.50  microns  < G )  and  0.85  microns  (U)  at 
2  0  deg  i  een  el*  vat  un  1 1.  t  he  solar  azimuth.  At  four  to  five 
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barely  recognizable.  Therefore,  the  intensity  decrease  at  G 
serves  as  stable  reference  for  that  at  U: 

LCR=log10  MU/G)4  degrees  /(U/G)0  degrees!. 

It  follows  that  for  periods  of  a  few  months  duration,  the 
upper  level  of  LCR  corresponds  to  twilights  undisturbed  by 
distant  clouds  (whether  volcanic  or  not),  while  reduced  LCR 
could  indicate  such  cloud  disturbances.  Because  the  neutral 
point  observations  were  made  in  white  light  (  \  ^  =  0.55 

microns),  it  is  unlikely  that  they  will  respond  to  the  cloud 
disturbance  indicated  by  the  LCR  values. 
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